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Abstract 

We report the first observation of the decays —)■ pADU')~. The data sample of 711 fb“^ 

used in this analysis corresponds to 772 million BB pairs, collected at the T(4S') resonance by 
the Belle detector at the KEKB asymmetric-energy e"’'e“ collider. We observe 19.8 cj and lO.Scj 
excesses of events for the two decay modes and measure the branching fractions of B^ —)■ pAD~ 
and B^ —)■ pAD*~ to be (25.1 ± 2.6 ± 3.5) x 10“® and (33.6 ± 6.3 ± 4.4) x 10“®, respectively, 
where the first uncertainties are statistical and the second are systematic. These results are not 
compatible with the predictions based on the generalized factorization approach. In addition, a 
threshold enhancement in the di-baryon (pA) system is observed, consistent with that observed in 
similar B decays. 

PACS numbers: 13.25.Hw, ll.30.Er, 12.15.Hh 
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In the years since the ARGUS and CLEO collaboration Erst observed baryonic B decays 
BEI, many three-body baryonic B decays {B —> BB'M) have been found |3H7j, where 
BB' denotes a baryon-antibaryon system and M stands for a meson. Although the general 
pattern of these decays can be understood as the interplay between the short-distance weak 
interaction and the long-distance strong interaction P], theories still have difficulties adjust¬ 
ing for various details such as the angular correlation between the energetic outgoing meson 
and one specihc baryon (B) in the di-baryon system HIEHII]. 

A popular theoretical approach used to investigate the three-body baryonic decays is 
generalized factorization. This method smears the correlation between the weak decay and 
the fragmentation and allows B — BB'Mc decays (with Me denoting a charmed meson) to 
be categorized into three types: current type, where the BB' pair is formed by an external 
W with other quarks; transition type, where the W is internal and forms BMc; and hybrid 
(current-|-transition) type [12]. The B^ pAD^*^~ [T2| decay belongs to the first type 
whereas its corresponding charged mode, B~^ —is of the last type. Using this 
approach. Ref. |T2| predicts the branching fractions 

i3(R° ^ pAD~) = (3.4 ± 0.2) x 10“®, 
i3(R® ^ pAD*-) = (11.9 ± 0.5) X 10“®, 

( 1 ) 

B{B+ pAD°) = (11.4 ± 2.6) x 10“®, 
i3(R+ ^ pAD*^) = (32.3 ± 3.2) x 10“®. 

There are two salient features of the predicted results. First, the ratios of the branching 
fractions of the decays into D* to the analogous decays into D are Ri 3 : 1. Secondly, the 
branching fraction of the hybrid-type decay is also ~ 3 times larger than the corresponding 
current-type decay. The measured branching fraction for B~^ —)■ pAD^ is consistent with the 
theoretical calculation based on the factorization approach [121 El. 

In most B BB'M decay studies, the hnal-state di-baryon system is observed to favor 
a mass near threshold EEHIZ]. While this “threshold enhancement effect” is intuitively 
understood in terms of the factorization approach, such enhancements are not seen in B~^ —)■ 
pAJ/ip nor R+ —)■ A+A“A'+ [lEl fTH] . More intriguingly, the factorization approach fails 
to provide a satisfactory explanation for the M-p angular correlations in B~ —)■ ppK~, 
R® —)• pAvr", and B~ —)■ ppD~ [71 [9IITT]. A striking difference between the non-zero angular 
asymmetries of R“ —)■ ppD*~ and R“ —)■ ppD~ was also reported in Ref. [Slini, for which 
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a theoretical explanation was attempted in Ref. [20] • A study of pure current-type decays 
like —)■ pAD^*^~ is useful to shed more light on the afore mentioned phenomena. In this 
paper, we report the hrst observation of —)■ pAD^*^~ decays using data from the Belle 
experiment. 

The data sample used in this study corresponds to an integrated luminosity of 711 fb“^ 
or 772 X 10® BB pairs produced at the T(4S') resonance. The Belle detector is located at 
the interaction point (IP) of the KEKB asymmetric-energy e~^ (3.5 GeV) e~ (8 GeV) col¬ 
lider [211122]. It is a large-solid-angle spectrometer comprising six specialized sub-detectors: 
the Silicon Vertex Detector (SVD), the 50-layer Gentral Drift Ghamber (GDG), the Aerogel 
Gherenkov Gounter (AGG), the Time-Of-Flight scintillation counter (TOP), the electro¬ 
magnetic calorimeter, and the and muon detector (KLM). A superconducting solenoid 
surrounding all but the KLM produces a 1.5 T magnetic held [2S1 124] . 

The hnal-state charged particles, tt^, and ^p\ are selected using the likelihood infor¬ 
mation from the combined tracking (SVD, GDG) and charged-hadron identihcation (GDG, 
AGG, TOP) systems [2S]. The R® —)• pARb)- signals are reconstructed through the sub¬ 
decays D~ —j- iP+TT^TT", D*~ —j- R® —)■ and A —)■ pTr~^. The distance of 

closest approach to the IP by each charged track is required to be less than 3.0 cm along the 
positron beam [z axis) and 0.3 cm in the transverse plane.The pion and kaon identihcation 
efficiencies are in the range of 85-95% while the probability of misidentifying one as the 
other is 10-20%, both depending on the momentum. The proton identihcation efficiency is 
90-95% for the typical momenta in this study, and the probability of misidentifying a proton 
as a pion (kaon) is less than 5% (10%). The candidate A is required to have a displaced 
vertex that is consistent with a long-lived particle originating from the IP and an invariant 
mass between 1.102 and 1.130 GeV/c^. The particle-identihcation criterion is omitted for 
the daughter pion in the A reconstruction due to the low background rate. For a R®, we 
require the reconstructed invariant mass to lie between 1.72 and 2.02 GeV/c^. For D~ and 
D*~, we require \Md- — 1870 MeV/c^| < 10 MeV/c^, \Md*- — 2010 MeV/c^| < 150 MeV/c^, 
and \Mf),- — Mp,o — 145 MeV/c^| < 9 MeV/c^, where and M^o are the reconstructed 

masses of R^*^“ and R®, respectively. 

We identify the signals using two kinematic variables: the energy diherence (AR) and 


6 


the beam-energy-constrained mass (Mbc), 


E^ -E'beam 

: - , P) 

Mbc = , 

where Eb and pb are the energy and momentum of the B meson and -Ebeam is the beam 
energy, all measured in the T(4S') center-of-mass (CM) frame. 

We optimize all selection criteria using Monte Carlo (MC) event samples before examining 
the data. These samples, both for signal and background, are generated using EvtGen [2S] 
and later processed with a GEANT3-based detector simulation program that provides the 
detector-level information j27j . 

Using the generated MG samples, the £t region is dehned as —0.1 GeV < AE < 0.3 GeV 
and 5.22 GeV/c^ < Mbc < 5.30 GeV/c^ while the signal region is dehned as \AE\ < 0.05 GeV 
and 5.27 GeV/c^ < Mbc < 5.29 GeV/cT 

Two major sources contribute as background: qq {q = u, d, s, c) production, 

also known as the continuum background, and other b ^ c dominated B meson decays, 
labeled generically as B decays in this paper. 

To suppress the continuum background, we use the difference between its jet-like topology 
and the spherical i?-decay topology. We calculate the distributions of 23 modihed Fox- 
Wolfram moments from the hnal-state particle momenta given by the signal and background 
MG [2H1 ESj- A Fisher discriminant that enhances the signal and background separation 
with a weighted linear combination of the moments is then calculated EDI. We augment 
the obtained probability density functions (PDFs) of the Fisher discriminant for the signal 
and background with two more variables to form the signal (background) likelihood £s(b): 
the axial distance (Az) between the vertices of the candidate B and the remaining hnal- 
state particles — presumably from the other B — and the cosine of the polar angle of the 
B momentum (cosOb) in the GM frame. The PDFs used for the modihed Fox-Wolfram 
moments, A^, and cos6b are bifurcated Gaussian functions, the sums of three Gaussian 
functions, and second-order polynomials, respectively. 

To suppress the background, we optimize the selection criteria for [£s/('Cs + 'Cb)]_d(d*) < 
aD(D*), \Mb- — 1870 MeV/c^l < ( 3 d MeV/c^, and \Md*- — Mqo — 145 MeV/c^| < ( 3 d* 
MeV/c^ simultaneously and obtain aD = 0.53, aD* = 0.40, ( 3 d = 10, and ( 3 d* = 9. The 
(3 selections correspond to ±2.4cr and ±12.4(j selections around the nominal Md*- and 
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Mjj*- — Mf) 0 . This procedure maximizes the hgure of merit, Ns/\/Ns + Nb , where Ns and 
A^b are the expected yields of signal and background, respectively, in the signal region. We 
use the theoretical expectations in Eq. (1) to obtain Ns and normalize the qq and generic 
B MC samples to the integrated luminosity to obtain A^b- After applying all the selection 
criteria, the fractions of events with multiple signal candidates are found to be 3.5% and 
5.6% in the D and D* modes, respectively. To ensure that no event has multiple entries in 
the fit region, we retain the B candidate with the smallest vertex-fit iri each event, where 
the vertex-£t is performed using all charged tracks from the B candidate except those from 
A. 

We model the signal AE distribution with the sum of three Gaussian functions; and the 
Mbc distribution with the sum of two Gaussian functions. We model the background AE 
shape with a second-order polynomial; and the Mbc shape with an ARGUS function 1311. 
We determine the PDF shapes with MG samples and calibrate the means and widths of the 
signal PDFs using a large control sample of B^ —decays from the data. The 
signal yields are extracted separately from eight di-baryon (pA) invariant mass bins, in the 
ranges of 2.05-3.41 GeV/c^ for the D mode and 2.05-3.30 GeV/c^ for the D* mode. We 
obtain the signal using a two-dimensional extended unbinned maximum likelihood fit in AE 
and Mbc. 

Figure [^illustrates the fit results of the lowest and highest pA mass bins for the D and 
D* modes. We observe clear signal peaks with very low background in the lowest bin, 
indicating an enhancement near threshold. As the efficiency is dependent on M^^, Table [^ 
lists the efficiencies and fitted yields in all mass bins for the two modes. Note that the 
efficiencies shown do not include the sub-decay branching fractions. 

Assuming that the branching fractions of T(4S') decaying to the charged and neutral BB 
pairs are equal, we use the efficiency and fitted yield in each mass bin to calculate the differ¬ 
ential branching fraction and integrate over the entire mass range to obtain the branching 
fraction B = Aj/ej)/(n-Bsubdecay x A^sP x Gpid), where i is the mass bin number, W and 
e* are the bin-dependent fitted yield and selection efficiency, respectively, iSsubdecay and 
are the sub-decay branching fraction and the number of BB pairs, respectively, and Gpid is 
the charged-particle identification efficiency correction between MG and data (0.92 for the D 
mode and 0.85 for the D* mode). Figure [^ shows the results, where both modes have visible 
peaks near threshold. The data are fit with an empirical threshold yield, x 




FIG. 1. Projections of typical AE-M]^^ fits to data for events in the signal region of the orthogonal 
variable. The peaking and flat red dotted lines represent the signal and background components; 
the blue solid lines with the dotted areas represent the combined PDFs with their la uncertainty 
bands. The top (bottom) four panels from top to bottom show the fits in the lowest and highest 
bin in the D (D*) mode. 


vs. the mass excess m = — Mp by varying a, b, c, and d. The obtained branching 
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FIG. 2. Differential branching fractions of the D (left) and D* (right) modes in Fit curves 

are based on an empirical threshold function (see text). 

fractions are: 

BiB^ ^ pKD-) = (25.1 ± 2.6 ± 3.5) x 10“®, 19.8 ct, 

(3) 

B{B^ ^ pKD*-) = (33.6 ± 6.3 ± 4.4) x 10“®, lO.Sa, 

where the quoted uncertainties are statistical and systematic (described later), respectively, 
and the signihcance is estimated by the Z-score of the p-value for = 2 yi,- ln(L ma^ ,/ Lm) 
with 8 or 6 degrees of freedom representing the number of bins. Lmax and Lq are the 
likelihood values with and without the signal component in the £t, respectively, and i is 

TABLE I. The fitted signal yield and efficiency in each bin. To obtain a stable fit, we combine 
the last three bins in the D* mode into the sixth bin. 


(GeV/c2) 

D mode 

(GeV/c2) 

D* mode 

Yield Eff.(%) 

Yield Eff.(%) 

2.05-2.22 

57 ±8 12.2 ±0.0 

2.05-2.21 

19 ±5 12.2 ±0.0 

2.22-2.39 

24 ±5 10.5 ±0.0 

2.21-2.36 

9 ±3 10.2 ±0.0 

2.39-2.56 

14 ±4 9.5 ±0.1 

2.36-2.52 

5 ±3 8.7 ±0.0 

2.56-2.73 

8±3 9.8 ±0.1 

2.52-2.68 

2±1 8.4 ±0.1 

2.73-2.90 

3±2 10.4 ±0.1 

2.68-2.83 

3 ±2 7.6 ±0.1 

2.90-3.07 

7±3 10.9 ±0.2 

2.83-3.30 

1±1 6.3 ±0.1 

3.07-3.24 

1±2 10.8 ±0.3 



3.24-3.41 

2 ±2 11.4 ±0.7 



Total 

117± 12 


39 ±7 
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again the mass bin index. The measured branching fractions are clearly incompatible with 
the theoretical predictions for both the D and D* modes [12]. This indicates that the model 
parameters used in the calculation need to be revised and, perhaps, some modification of 
the theoretical framework is required. 

To extract the decay angular distributions, we divide cos9p^(*) into eight bins, where 
6pjj(*) is dehned as the angle between the proton and meson directions in the pA rest frame. 
We follow the same procedure to determine the differential branching fractions in cos6pj^(*) 
as in determining those in Table [IT| lists the htted signal yields and efficiencies in the 

cosOpDi*) bins; Fig. shows the differential branching fractions. The efficiency is determined 
with the MC sample, including the threshold enhancement effect as observed in the data. 

We define the angular asymmetry Aq = , where represents the branching 

fraction of positive (negative) cosine value. The results are 


Ag{B^ pAD-) = -0.08 ± 0.10, 
Ae{B^ pAD*-) = +0.55 ± 0.17, 


(4) 


where the uncertainty is purely statistical since the correlated systematic uncertainties cancel 
in the Ag calculation. The angular distributions of the D and D* modes appear to have 
distinct trends, even though they are both categorized as current-type decays. More data 
are needed to make the result conclusive. 


TABLE II. The fitted signal yield and efficiency in each cos6pjj(,) bin. 


cosOpj^M 

D mode 

D* 

mode 

Yield 

EfF.(%) 

Yield Eff.(%) 

-1.00 

-0.75 

10 + 4 

9.0 

3 + 2 

8.6 

-0.75 

-0.50 

17 + 5 

10.5 

1 ± 1 

10.2 

-0.50 

-0.25 

16 + 4 

11.5 

1 ± 1 

11.3 

-0.25 

-0.00 

15 + 4 

12.2 

2 + 2 

12.2 

+0.00 

- +0.25 

19 + 5 

12.8 

7 + 3 

12.7 

+0.25 

- +0.50 

15 + 4 

13.0 

7 + 3 

13.0 

+0.50 

- +0.75 

16 + 5 

12.6 

9 + 3 

12.8 

+0.75 

- +1.00 

7 + 3 

11.5 

8 + 3 

11.5 
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FIG. 3. Differential branching fractions of the D (left) and D* (right) modes in cos9p^(t,). The fit 
curves are second-order polynomials, as suggested by Ref. [20] . 

Three major categories of systematic uncertainties are considered: in the signal yield de¬ 
termination, in the efficiency estimation, and in translating the signal yields and efficiencies 
into the branching fractions. Table |III| lists all the systematic uncertainties. 

We observe a mild peaking background in the Mbc fit region due to —)■ pAD*^, 

plausibly by the replacement of the low-momentum 7r° in D*^ —with an unaffiliated 
7r“ or K~ to reconstruct a D*~. To study its contribution to the uncertainty in the D* 
mode, a dedicated MC sample of this background mode is generated. Based on its current 
branching fraction upper limit [TT], we subtract 0.5 events from the extracted signal yield and 
assign ±0.5 events as the systematic uncertainty. We have verihed that our signal extraction 
method is robust and see negligible systematic bias in the signal yield when assuming 0.1 to 
10 times the theoretical branching fractions (about 1.6 to 160 events) in an MC ensemble 
test. 

For the reconstruction efficiency, we consider the following systematic uncertainties: the 
signal MC modeling for the threshold enhancement effect using the bound state assumption, 
charged track reconstruction, charged hadron identihcation, A reconstruction, background 
discrimination selections, and the PDF shapes. The modeling uncertainty is estimated 
by comparing the efficiency calculation based on two different MC samples, one generated 
assuming p-A bound states and the other with three-body phase-space decays, in each 
bin. As the result is highly threshold-enhanced, we use the efficiency given by the bound- 
state model to calculate the branching fractions and take the differences as the systematic 
uncertainties between the two models. The uncertainty is about 3 (2)% in the D {D*) 
mode, depending on the bins. For each charged track except the low-momentum pion in 
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D*~ —)• D^7r~, a 0.35% uncertainty is assigned to take into account the data-MC difference 
in the charged track reconstruction. For the low-momentum pion, a 2.5% uncertainty is 
assigned. We use the A —)• pn~ and D*~^ —)• —)■ K~tt^ samples to calibrate the 

MC TT^ identihcation efficiencies and assign uncertainties. For the A reconstruction, 

we estimate the uncertainty by considering the data-MC difference of tracks displaced from 
the IP, the A proper time, and A mass distributions. The uncertainties due to the «£,(*) 
selections are estimated separately with the control sample mode, —)• . We 

compare the data-MC efficiency differences with or without the a selections, where the non- 
negligible statistical uncertainties are also included. In both cases, the obtained —)• 

is found to be consistent with the world average, indicating overall reliability 
of our methodology. For the jSo and Pd* selections, we compare the widths of the peaking 
components in M^- and M^*- — Mjjo in the MC and data and quote the differences as the 
uncertainties. We also relax the shape variables of the signal PDF when htting the control 
sample and compare the difference to MC-determined PDF. The resulting difference in the 
calculated B{B'^ —)• KgD^*'>~) is negligible. 

TABLE III. The systematic uncertainties in the D and D* modes. The ~ signs indicate the 
dependence of the uncertainty. 


Item 

Systematic uncertainty (%) 

D mode 

D* mode 

Yield bias 

negligible 

1.3 (0.5 evt.) 

Modeling 

3 

2 

Charged track 

2.1 

4.3 

Charged hadron identification 

1.3 

1.8 

A identification 

4.0 

4.4 

Md-, Md*-—M^o window 

2.0 

negligible 

Ts/(Ts + Tb) requirement 

11.5 

11.0 

PDF shape 

negligible 

negligible 

^BB 

1.4 

1.4 

Sub-decay B 

2.2 

1.7 

Overall 

13.9 

13.1 
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In the translation from signal yields to branching fractions, we consider the uncertainties 
of -Bsubdecay and N^j^. The uncertainties of iSsubdecay are obtained from Ref. [3]. For Nsb, 
on- and off- resonance di-lepton events, —)■ qq MC and data difference, primary vertex 

sideband data, and statistical uncertainty are combined to estimate the uncertainty. 

In this paper, we have reported the hrst observation of the —)■ pAD~ and —)■ pAD*~ 
decays with branching fractions (25.1 ± 2.6 ± 3.5) x 10“® (19.8cr) and (33.6 ± 6.3 ± 4.4) x 
10“® (lO.Scr). The threshold enhancement effect observed in is found to be consistent 
with many other three-body baryonic B decays. The obtained branching fractions disagree 
with predictions based on the factorization approach, as do the measured ratios of branching 
fractions, both for the D and D* modes and for the charged and neutral B modes. We also 
hnd potential angular asymmetry in the D* mode but not in the D mode. Theoretical 
explanations, as well as conhrmation from experiments with sizable data sets, such as LHCb 
and Belle II, will be needed in the future. 
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